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ABSTRACT 

Aii — ice-ocean  feedback  mechanisms,  which  are  not 
conventionally  incorporated  within  either  climate  or  glacial 
models,  are  investigated  to  illustrate  their  potential  role  in 
generating  ice  advance/retrcat  on  the  time  scale  of  1 03— 104 
years;  i.e.  for  examining  the  internal  causes  for  the  ice 
oscillation. 

Three  main  feedback  loops  are  found  from  a  coupled 
air-ice-ocean  model  developed  in  this  paper;  (a)  ice  advance 
-  lower  air  temperature  -*  ice  freezing  -  ice  advance:  and 
(b)  ice  advance  -  higher  ocean  temperature  -  ice  melting  - 
ice  retreat;  (c)  ice  advance/retreat  -  modification  of  evap¬ 
oration  rate  -  change  of  ice  accumulation  rate  and  sea-level 
height  -  ice  advance/retreat  The  relative  strength  of  the 
three  feedback  mechanisms  determines  the  characteristics  of 
the  modes:  growing  or  decaying,  oscillatory  or  non- 
oscillatory.  The  solutions  show  the  generation  of  growing 
oscillatory  modes  with  the  time  scale  of  I03-104  years  in 
certain  parameter  ranges. 
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Fig.  I.  Physical  processes  in  the  coupled  air-ice-occan 
system. 


INTRODUCTION 


MODE!.  DESCRIPTION 


In  trying  to  understand  ice  oscillation  and  related 
climatic  change,  one  may  First  ask:  is  ice  advance/retreat 
caused  by  external  or  internal  forcings?  The  externa! 
forcings  are  agencies  outside  the  Earth  system,  such  as 
astronomically  caused  changes  in  the  Earth's  orbital 
parameters. 

Internal  forcings,  on  the  other  h3nd.  are  processes  and 
mechanisms  present  within  the  Earth  system,  giving  rise  to 
self-excited  glacial  fluctuations.  The  temperature-albedo 
feedback  mechanism  associated  with  long  and  short  wave 
radiation  fluxes  at  the  Earth's  surface  are  given  considerable 
attention  in  glacial  and  climate  models  (Kailcnnnd  others. 
1979;  Oerlemans  and  van  der  Veen.  1984).  However,  the 
role  of  the  hydrological  cycle  on  the  ice  advance/retrcat  is 
given  less  attention  although  it  is  realized  that  ice 
freeze, 'melt  and  evaporation/precipitation  are  contributors  to 
the  change  Of  air  and  ocean  temperatures,  which  in  turn 
influence  the  ice.  The  purpose  of  this  paper  is  to 
investigate,  in  its  simplest  form,  a  possible  feedback 
mechanism  in  the  air-ice-ocean  system  regarding  the 
hydrological  cycle.  The  ice  advance/  retreat  on  the  continent 
is  not  included  in  this  work. 

Figure  I  shows  the  main  physical  processes  in  the 
coupled  system.  From  the  global  point  of  view,  ice  advance 
isolates  the  ocean  from  the  atmosphere  and  reduces  the  heat 
loss  from  the  ocean  surface,  which  in  turn  decreases  ihe  air 
temperature  and  increases  the  ocean  temperature.  The  air 
and  ocean  feed  back  into  the  ice  through  two  different 
mechanisms  (Fig.  2):  (I)  thermodynamical  feedback;  the 
increase  of  ocean  temperature  melts  the  ice  and  causes  ice 
retreat  (negative  feedback  between  ice  and  ocean),  however, 
the  decrease  of  air  temperature  causes  freezing  and  makes 
the  ice  further  advance  (positive  feedback  between  ice  and 
air):  (2)  mechanical  feedback;  the  modification  of  the 
atmospheric  and  oceanic  temperatures  varies  the  surface 
evaporation  rate,  which  in  turn  changes  the  se3-level  height 
and  the  ice  accumulation  rate.  Both  effects  lead  to  a 
further  change  of  the  ice  flow. 


Ice  advance/ retreat 

The  total  surface,  ice  covered,  and  open  ocean  areas 
are  set  to  be  1.  <Vj,  and  1  -  A’j,  respectively  (Fig.  I).  The 
ice  is  assumed  to  be  isotropic  with  thickness  Hy  Following 
Weertman  (1957).  the  ice-spreading  rate  is  computed  by: 

<ih>  =  4^-pig//,]"  (I) 

where  p;  is  the  ice  density.  Hx  is  the  height  between  the 
ice  top  and  the  sea-level  (Fig.  I).  Values  of  n  vary  from 
about  1.5  to  4.2  with  a  mean  of  about  3.  and  for  randomly 
oriented  polycrystalline  ice  at  -10  C  and  n  -  3,  a  value 
A  =  3  x  10  *  a'1  kPa  3  is  reasonable  for  ice  stress 
(Weertman,  1973).  The  angle  brackets  denote  values  averaged 
over  ice  thickness.  The  time  rate  of  change  of  ice  coverage 
.Vj  contains  mechanical  and  thermodynamical  (ice  edge 
freezing/melting)  parts: 


Fig.  2.  Feedback  mechanisms  among  air,  ice.  and  ocean 
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where  //;  is  ihe  ice  thickness,  and  F  is  the  freezing  rate  at 
the  ice  edge.  Following  Killworth  (1979),  the  heat  flux  is 
accomplished  through  an  exchange  coefficient  K  and  is  of 
the  form  K(T  -  Tf).  where  Tf  represents  the  freezing  point 
of  sea  water.  Therefore  the  freezing  rate  is  computed  by: 


F  =  ~  7y~  iPwcpwAwi(rw  -  rf)  +  PncpafCai(Ta  -  Tf)l 
H  -l 

(3) 

where  Kwr  A'aj  arc  the  heat  exchange  coefficients  for 
water-ice  and  air-ice.  respectively.  L-t  is  the  talent  heat  of 
ice.  If  Tw.  Ta  are  greater  than  Tf,  F  is  taken  as  a 
negative  value,  which  means  ice  melting. 

Time  rale  of  change  of  ice  thickness 

Based  on  the  equation  of  continuity.  Shumskiy  (1965) 
proposed  a  method  of  determining  ice-thickness  changes  by 
comparing  the  accumulation  rate  with  horizontal  strain-rates. 
The  time  rate  of  change  of  ice  thickness  becomes: 

3"i 

—  =  CS  -  -  Fj  W; 


where  cs  is  ice  accumulation  rate,  and  T;  is  vertical  mean 
ice  velocity,  which  links  to  ice  thickness  and  bedrock 
topography.  It  is  commonly  believed  that  the  non-linear 
inertial  term  — Tj-v//j  in  ice-sheet  models  relates  to  the  long 
time  scale  of  glacial  cycles;  however,  for  simplicity  this 
term  is  neglected  from  the  present  model,  i.e.: 


f^i 

d: 


//i<ch> 


(4) 


which  means  that  the  current  coupled  system  excludes  the 
variations  with  long  time  scales  comparable  to  the  glacial 
cycles. 

Time  rate  of  change  of  ocean  depth 

From  a  thermodynamical  point  of  view,  the  surface 
evaporation  and  ice  freezing/melting  are  two  major  processes 
to  change  the  ocean  elevation,  i.e.: 

3HW  (I  -  Ufa F 


where  Hy.  is  the  depth  of  the  ocean.  pw  is  the  sea-water 
density,  and  E  is  the  surface  evaporation  rate  computed  by 
the  bulk  formulx 


*  -  <»  -  ^Dl^alr^wJ  -  *ai 

P  w 

where  CD  is  the  drag  coefficient.  I'a  is  the  surface  wind, 
is  the  saturated  mixing  ratio,  and  <?a  is  mixing  ratio  near 
the  ocean  surface. 

Heat  equations  for  air  and  ocean 

The  heat  equations  for  the  atmosphere  and  ocean  are: 


MEAN  STATE  AND  PERTURBATIONS 

Paramctcrigulion  of  perturbations 

Mean  state  of  the  coupled  system  (Ay  H-r  ffw,  7*a, 
7*w,  E)  is  obtained  by  the  zero  time  rate  of  change  in  the 
basic  Equations  (2),  (4),  (5),  (6)  and  (7).  When  the  coupled 
system  is  perturbed  from  its  mean  state,  the  feedback 
mechanisms  make  the  perturbation  either  grow  or  dampen. 
Assume  that  the  perturbation  of  the  ocean-surface 
evaporation  rate  is  caused  by  the  fluctuations  of  sea-surface 
temperature  (7"  w),  surface  wind  speed  (V  a),  and  the  ice 
coverage  (A’' ;,),  i.e.: 


Pad9s(7"w) 

£■=(.-  Ayc^l-^-Ty  ♦ 

(1) 


E 

- v  i  • 

(1  ~  A’;)  1 


(3) 


(8a) 


Ice  advance  results  in  an  increase  of  the  ocean  temperature 
and  a  decrease  of  the  surface  air  temperature.  The  greater 
ice  coverage  will  partly  insulate  the  atmosphere  from  the 
ocean,  decreasing  the  evaporation  (see  term  (3));  but  in 
locations  where  the  ocean  and  atmosphere  are  in  direct 
contact,  the  warmer  ocean  temperatures  will  increase  the 
ocean  evaporation  rate,  which  is  called  the  local 
SST-evaporation  coupling  (see  term  (1)).  Also  the  colder  air 
temperatures  will  strengthen  the  thermally  driven  surface 
wind  (due  to  the  larger  surface  air  temperature  gradient) 
and  consequently  increase  the  evaporation  rate,  which  is 
called  the  wind -evaporation  coupling  (see  term  (2)).  The  net 
effect  as  to  whether  the  ice  advance  increases  or  decreases 
the  evaporation  rate  is  thus  uncertain.  According  to  Chu's 
(1989)  thermally  driven  surface-wind  formula,  the  term  (2) 
in  Equation  (8a)  is  an  order  of  magnitude  smaller  than  the 
term  (1),  and  therefore  is  neglected  from  the  equation, 
leading  to: 


E'  =  (I  -  <V;)7 


Pa  dgs<rw> 
Pw  <ir 


-T 
‘  w 


(8b) 


where  y  (=  CD  |  t'a  j )  is  the  moisture  exchange  coefficient. 

The  relationships  among  the  surface  evaporation  rate 
£' ,  the  precipitation  rate  on  the  ocean  Aj..  and  the  ice 
accumulation  rate  cj  are  simplified  as: 

?\  =  (i  -  Ayr,  -  Ayr  (9) 

which  means  the  surface  evaporation  is  the  source  for 
precipitation  and  ice  accumulation.  Following  Killworth 
(1979),  the  heat  fluxes  are  computed  by: 
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where  Cpa.  cpw  are  the  specific  heats  for  air  and  ocean, 
respectively,  and  Aaj,  A‘wa.  and  A.'w  -  are  heat-exchange 
coefficients  between  the  corresponding  two  components. 

The  perturbation  of  lateral  freezing  rate  is  computed 

by: 


where  7\.  Tw  are  temperatures  of  air  and  ocean.  7/a  is  the 
depth  of  the  atmosphere.  <7a],  Qw j  are,  respectively,  the 
heat  fluxes  from  air  and  ocean  to  ice  and  <?wa  is  the  heat 
flux  from  ocean  to  atmosphere.  Equations  (2),  (4),  (5),  (6), 
and  (7)  are  the  basic  equations  of  our  model. 


(Pw^pw'A  wi^w  +  P3cpaAai7V 

pj^i 
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where  £j  is  the  latent  heat  of  the  ice. 
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Timescales  of  perturbations 

Four  main  timescales  appear  in  this  coupled  system. 
They  are: 


*ai  *wa 
- ,Vj  +  - (I  —  ft’;) 

L"a  '  V 


(12) 


which  is  the  atmospheric  relaxation  timescale: 

A  wi  PaCp3^W3 

r„,  =  - .V:  +  - (I  -  .V:) 

„  w  pw  w 

which  is  the  oceanic  relaxation  timescale. 

<?wa  +  0  "  AfjMT^lr* 

7e  S - - - 

■@wa  Tw^wa2r  Pa  3?s(fw)  l 


(13) 


Pwcpw77w  Fw  ^  '  Afj) 


(14) 

which  is  the  timescale  for  the  change  of  evaporation  rate 
due  to  the  ocean  warming/cooling  induced  by  the  ice 
advance/retreat,  and: 


«  k^jt  ns) 

which  is  the  ice-spreading  timescale. 

By  using  the  values  listed  in  Table  I,  these  four 
timescales  are  estimated  as: 


ra  -  1 1 .6  d,  rw  -  66  year.  Ts  -  3000  year 
(for  /7,  =  10  m),  rc  -  4660  year 
where  we  assume  that  |£?wa|»|(?wj|- 


(16) 


TABLE  I.  MODEL  PARAMETERS 


*ai  *  3  x 

I0's  m  s'1 

Kwa  =  10‘*  m  s'1 

*wi  “  4  * 

10"5  m  s1 

y  =  1.5  x  Kwa  (Gill.  1982) 

Uosbcrger, 

1987) 

f/2  =  8  km 

Hw  =  5  km 

f/j  =  2  km 
/V:  =  0.2 

(?wa  (Parpa/;(j*  =  1  *5  s  .  , 

«  =  3,  .4  =  j  *  10’8  year'1  kPa'* 

I  *v]*i  on 

9/lw 

—  =  1  -  (1  -  .V;)xor]Arj  (21) 

where  x  s  pj//j /(pw//w),  and  o  is  a  non-dimensional 
parameter  defined  by: 

a  -  7  w7  w  —  Og/lgTg)  .  (-2) 

For  a  >  0,  the  ice  has  a  "restoring  force’  induced  by  the 
ice  advance/retreat.  The  mechanical  feedback  of  the  air  and 
ocean  into  the  ice  can  be  easily  noticed  from  the  feedback 
of  hj  and  h'w  change  due  to  ft'j.  into  the  ice 

advance/ retreat  through  the  first  term  in  the  righthand  side 
of  Equation  (19).  This  feedback  largely  depends  on  the 

parameter  v  s  Hw/Hy  The  larger  the  value  of  v.  the 

stronger  this  feedback  mechanism  will  be.  Figure  3  shows 

the  relationships  among  the  model  variables. 


Fig.  3  Relationships  among  the  model  variables. 


The  general  solutions  of  Equations  (19).  (20)  and  (21) 
have  the  following  forms: 


Afjfi)  =  ZbjC  J  . 


h-U) 


_  W / 

I Xj-e  1  . 


h'(i ) 


Ujt  1 


(23) 


where  bj.  Cj.  and  dj  ( j  =  1.2,3)  are  the  integral  constants, 
and  w,,  u5.  and  u3  are  the  eigenvalues  which  are  the  roots 
of  the  third-order  algebraic  equation: 

us  +  o,u!  +  tfjW  +  c3  =  0 .  (24) 


Since  the  atmospheric  and  oceanic  relaxation  timescales 
Ta  and  7W  are  so  much  shorter  than  the  other  two 
timescales  7$  and  7e.  the  atmospheric  and  oceanic 
temperature  perturbations  Tj  and  T'w  almost  "instant¬ 
aneously"  follow  the  ice  advance/retreat  processes,  i.e.: 

7a  “  -v,aTaiVi-  rw  »  4w7w,Vj.  (17) 

where  .4a  =  «?wa  +  Oai)/(paCp3f/a).  =  Q wa/(Pwrp w^w )- 
Equation  (17)  clearly  shows  that  the  ice  advance 
(iVj  >  0)  causes  the  cooling  of  the  atmosphere  (7^  <  0)  and 
the  warming  of  the  ocean  ( T ^  >  0).  Substitution  of 
Equation  (17)  into  Equation  (111  leads  to  the  connection  of 
lateral  freezing  rate  F’  to  ft'j: 

F  =  -H-l(awAyMTw  —  (18) 

where  otg  =  paCpaA.'ij/(pj£|Wj),  =  pwCpWA.w/(p;E|//j) . 

Stability  parameters  a  and  v 

Linearization  of  the  basic  Equations  (2),  (4),  and  (5) 
leads  to: 

3A’|  (h[  -  k-w) 

IT  =  W* - irx - (l  ■  Vi)TeloVi  <19> 


RESULTS 

We  compute  all  roots  of  Equation  (24)  for  different 
values  of  the  parameters  o  and  v,  and  obtain  three  roots  at 
each  pair  of  the  parameters  (o,  v).  One  root  among  the 


e.o-. 


I - ! - 1 - 1 - i - 1 
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a 

Fig.  4.  Dependence  of  growth  rate  ur  on  stability 
parameters  a  and  v. 
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three  Has  a  negative  re3l  part  throughout  the  whole 
parameter  space,  representing  the  damping  modes,  in  which 
we  are  not  interested  here.  The  other  two  roots,  u,  and  u„. 
have  positive  real  parts  somewhere  in  the  parameter  space, 
representing  the  existence  of  growing  modes  in  certain 
parameter  ranges.  Figure  4  shows  the  dependence  of  growth 
rate  (10'15s*)  on  o  for  three  different  values  of  v.  The 
horizontal  axis  a  changes  from  0  to  10,  representing  the 
increase  of  thermodynamic  stability.  Oscillation  takes  place 
only  in  certain  range  of  a  {Table  11).  i.e.  for  oH  »  a  >  oL, 
Wj  *  0.  Out  of  this  range,  i.e.  for  a  <  crL,  or  o  >  oH. 
(J;  =  0,  which  represents  non-oscillatorv  modes.  Figure  5 
shows  the  dependence  of  the  period  (2n/Wj)  on  a  for  three 
different  values  of  v. 

'•21 


n 


Fig.  5.  Dependence  of  period  of  oscillatory  modes.  2ff/|Uj|, 
on  stability  parameters  a  and  v. 


TABLE  II.  THE  OSCILLATION  INTERVALS 


V 

°L 

aH 

50 

3.75 

5.75 

100 

3.75 

5.00 

200 

4.00 

4.75 

The  results  are  listed  below: 

(1)  The  two  parameters  cr  and  v  play  different  roles  in 
the  coupled  system.  The  negative  feedback  measured  by  a 
stabilizes  the  system;  however,  the  positive  feedback 
measured  by  v  destabilizes  the  system.  Because  of  the 
existence  of  positive  feedback  mechanisms,  the  coupled 
system  is  unstable  unless  the  parameter  a  exceeds  some 
critical  value  ac.  Figure  4  indicates  that  oc  =  4.25. 

(2)  Four  different  types  of  modes  are  generated  due  to 
the  relative  strength  of  positive  to  negative  mechanism: 

o  <  aL,  non-oscillatory  growing, 

i  o  i  ac,  oscillatory  growing, 

ac  S  a  i  cr,,.  oscillatory  damping, 

a  >  o(i.  non-oscillatory  damping. 


(3)  The  parameter  v  serves  as  an  instability  parameter. 
In  the  growing  regime  (Fig.  4).  the  larger  the  value  of  v. 
the  larger  the  growth  rate  will  be.  The  mean  growth  rate 
of  the  growing  modes  varies  from  2-6  x  lO'IJs'3.  The 
doubling  time  during  which  the  perturbation  doubles  its 
sirength  is: 


T a  =  -  ~  3660  -  1 1  000  year. 

ur 

(4)  The  parameter  v  3lso  largely  impacts  on  the 
periodicity  of  the  coupied  system.  First,  v  affects  the 
critical  values  of  oL  and  <rH.  The  larger  the  value  of  v.  the 
smaller  the  range  of  oscillatory  mode  will  be.  The  mean 
period  varies  f rom  2000-3000  years  (for  v  =  200)  to 
4000-12  000  years  (for  v  =  100). 


CONCLUSION 

This  conceptual  air-ice-ocean  coupled  model  shows  a 
possible  positive/negative  feed-back  mechanism,  induced  by 
the  hydrological  cycle,  among  the  atmosphere,  ice,  and 
ocean.  The  theory  predicts  the  generation  of  oscillations  in 
ice  coverage,  ice  thickness,  air  and  ocean  temperatures,  and 
sea-level  height,  on  the  time  scale  of  1 03 — 1  O'*  years. 
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